Intersubband-induced spin-orbit interaction in quantum wells 
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Recently, we have found an additional spin-orbit (SO) interaction in quantum wells with two subbands 
[Bernardes et ai, Phys. Rev. Lett. 99, 076603 (2007)]. This new SO term is non-zero even in symmet- 
ric geometries, as it arises from the intersubband coupling between confined states of distinct parities, and its 
strength is comparable to that of the ordinary Rashba. Starting from the 8x8 Kane model, here we present 
a detailed derivation of this new SO Hamiltonian and the corresponding SO coupling. In addition, within the 
self-consistent Hartree approximation, we calculate the strength of this new SO coupling for realistic symmetric 
modulation-doped wells with two subbands. We consider gated structures with either a constant areal electron 
density or a constant chemical potential. In the parameter range studied, both models give similar results. By 
considering the effects of an external applied bias, which breaks the structural inversion symmetry of the wells, 
we also calculate the strength of the resulting induced Rashba couplings within each subband. Interestingly, 
we find that for double wells the Rashba couplings for the first and second subbands interchange signs abruptly 
across the zero bias, while the intersubband SO coupling exhibits a resonant behavior near this symmetric 
configuration. For completeness we also determine the strength of the Dresselhaus couplings and find them 
essentially constant as function of the applied bias. 



1. INTRODUCTION 

The coupling between spatial and spin degrees of freedom 
in semiconductors provides an interesting possibility for co- 
herently manipulating the electron spin via its orbital (charge) 
motion. For instance, the proposal of Datta and DasU for a spin 
field-effect transistor highlights the use of the spin-orbitlSO) 
interaction of Rashba,LWU which is electrically tunableJJU to 
control - via spin rotation - the flow of electrons between fer- 
romagnetic source and drain. 

In addition to the Rashba SO coupling present in het- 
erostructures with structural inversion asymmetry in the con=. 
fining potential, there is the Dresselhaus SO interactionU 
present in both bulk and confined structures with inversion 
asymmetry in the underlying crystal lattice. These spin or- 
bit interactions have played an important role in the exciting 
field of semiconductor spintronics as they underlie a num- 
ber of interesting physical phenomena and potential spin- 
tronic applicationsLUU For instance, the effective zitterbe- 
wegung of spin-polarized wave packets injected into SO cou- 
pled two-dimensional (2D) electron gases is a very interest- 
ing possibility ! I I The interplay of the Rashba and Dres- 
selhaus interactions can give rise to conserved spin-rotation 
symmetriesUU relevant for devising-jpbust SO-based devices 
operating in the nonballistic regimeLJ 

Recently, a new type of SO interaction arising in quantum 
confined systems with two subbands has been foundLJ. Un- 
like the usual Rashba SO, this new SO term is nonzero even 
in wells with full structural inversion symmetry (and hence 
it does not produce spin splitting). This essentially follows 
from the distinct parities of the confined states (even and odd), 
which can couple via the derivative of a symmetric poten- 
tial. This intersubband-induced SO coupling is quadratic in 
the crystal momentum, unlike the Rashba and the (linearized) 
Dresselhaus terms in wells.LI As shown in Ref. \~j, this SO 
coupling can give rise to an unusual zitterbewegung (both in 
position and in spinl—l) and a nonzero spin Hall conductivity. 



Here we complement and extend the work of Ref. fL5|: (i) 
We present a more thorough derivation of the intersubband- 
induced,. SO interaction, starting from the 8x8 Kane 
mode]|_H_LJ within the k • p approach. We also slightly gen- 
eralize the derivation for confined systems with more than 
two subbands and structurally asymmetric potentials in which 
the usual Rashba-type SO interaction is present, (ii) We per- 
form a detailed investigation of the relevant SO couplings via 
a self-consistent scheme where we solve both Poisson and 
Schrodinger equations numerically (Numerov method) within 
the Hartree approximation. We consider realistic modulation- 
doped single and double quantum wells with applied external 
biases, which can change the spatial symmetry of the wells, 
and having either a constant areal electron density or a con- 
stant chemical potential. 

Our simulations focus on wells with two subbands. 
For nonzero applied biases we calculate not only the 
intersubband-induced SO coupling r\ but also the Rashba-type 
couplings ao, a\ for the first and second subbands, respec- 
tively. For completeness, we also calculata-the linearized 
Dresselhaus SO couplings for each subband. U For both the 
constant density and constant chemical-potential models con- 
sidered, we find sizable values of the intersubband SO cou- 
pling T] as compared to the usual Rashba and Dresselhaus 
couplings. Interestingly, for double wells near the symmetric 
(zero-bias) configuration we find that n has a resonant behav- 
ior, changing its magnitude by a factor of 10. On the other 
hand, the Rashba couplings for the first and second subbands 
abruptly change signs around the zero-bias voltage. The Dres- 
selhaus couplings do not show any noticeable behavior around 
this point, being essentially constant as a function of the ap- 
plied bias. 

We note that the SO coupling constants 77, ao, and a\ con- 
tain contributions from the potential well (and barrier) offsets, 
the electronic Hartree potential, and the external gate plus the 
modulation doping potentials. For the single wells investi- 
gated here, the external gate (+ modulation doping) is the 



dominant contribution to ao and a\, while r\ is mostly de- 
termined by the Hartree and structural offset contributions. 
On the other hand, for the double wells studied the electronic 
Hartree potential is the dominant contribution to 77, ao, and 
a\. Interestingly, the Hartree potential in this case is highly 
influenced by the external gate, particularly around the sym- 
metric (zero-bias) configuration, as the electrons can localize 
in either well for small (positive or negative) changes in the 
gate potential. This renders r\, ao, and ci more amenable to 
gate modulations in double wells as compared to single wells. 
Next we outline our work. 

In Sec. g we review the k • p approach and the Kane 
model. In Sec. [5] we present a detailed derivation of our effec- 
tive Hamiltonian for electrons in heterostructures with many 
confined states within the Kane model. In this section we 
also show the relevant expressions for the new intersubband- 
induced SO coupling r\ and those for the Rashba a (and Dres- 
selhaus /?) SO couplings as well. In Sec. we describe the 
quantum wells investigated and (briefly) the standard self- 
consistent calculation performed. We present our results and 
discussions in Sec. g. In this section we focus specifically on 
realistic single and double-well systems. Section H summa- 
rizes our work. In Appendices ^ and [b] we show details of 
our self-consistent scheme to solve the relevant Schrodinger 
and Poisson equations. 



2. k p APPROACH AND KANE MODEL 

Here we briefly review the kp approach and use it to obtain 
the 8x8 Kane model relevant for pur, derivation of the new 

!9fZ 



intersubband-induced SO coupling 



As usual, to solve Eq. (g), we expand u„k in terms of the 
eigenstates m/o at k — [i.e., W(k - 0) = 0] obtained from 



H(k = 0) uio(r) = s m u m (r), 



(5) 



where / = 1,2, . . .N (in principle, N — > 00) indexes the dis- 
crete set of levels at k = [note that Eq. (Q) contains the SO 
interaction, even though W(0) = 0]. Substituting 



u„k(r) = 2_j a ni{k) u m (r), 



(6) 



i=\ 



into Eq. (g) and projecting the resulting expression onto the 
ui>o(r) eigenstate, we fincO 



2^ [{en - e„k + ^— - )8 W 



2niQ ' 



+ <r|— * • p + -^-5* • tr x VV(r)|/)l a nl (k) = 0. (7) 



mo Am^c 1 



Here we use the notation (r\l) = um{r) and define 



'|A|/)= I d 3 ru; Au m , 



(I 



with A denoting a Hermitian operator. 



(8) 



2.1. Basics of the k ■ p method 



The single-particle Hamiltonian, [for an electron with mo- 
mentum p in a periodic potentialLJU V(r) with SO is 



H=^-+ V(r) + -^—o- x VV(r) • p, 
2mo Ami.c 1 



(1) 



where mo is the bare electron mass and <x is a vector oper- 
ator defined in terms of the Pauli matrices. With the help of 
Bloch's theorem 4i„k{r) - exp(/£ ■ r) M„t(r) [M„t(r) has the pe- 
riodicity of the underlying Bravais lattice] we can rewrite the 
Schrodinger's equation Hifr„k = e n k4>nk, where n indexes the 
distinct solutions for each k vector, in the form 



[H(k = 0) + W(k)]u nk (r) = [e nk 



H 2 k 2 
2mo 



u„ k (r), (2) 



with 

H(k = 0) 



V 2 + V(r) + -^-<rxVV(r)./,, (3) 

2m Q 4mf.c z 



W(k) = —k ■ [p + -a x VV(r) 

mo \ 4moc z 



(4) 



2.2. 8x8 Kane model - bulk case 

As usual, in order to solve Eq. (g) we have to truncate the 
basis set by considering a finite number Af of zone-center ba- 
sis functions uio(r). In addition, since the k - Hamiltonian 
[Eq. (B|)] contains a SO term, it is convenient to choose linear 
combination of basis functions which are eigenstates of the to- 
tal angular momentum J = L + S, and its z component J z ; here 
L and S denote the orbital and spin angular momenta, respec- 
tively. In II- VI and III-V (both zincblend) compounds the rel- 
evant conduction and valence bands arise from the "bonding" 
p-type and "anti-bonding" s-type states, respectively. Follow- 
ing the notation of Refs. rj and \~\ we summarize in Table | 
the set of eight zone-center wave functions we consider here 
(the kets \JJ Z ) are also shown), which are the eigenstates of the 
zone-center Schrodinger's equation (|5|) for the periodic part of 
the Bloch function. Note that we use the standard state vector 
notation \S), \X), \Y), and |Z) to denote the symmetry of the 
corresponding "atomic orbitals" (tight-binding view). 

Using the ordered basis states U\, . . . , u% in Table ffl we can 
easily write out the matrix Hamiltonian [Eq. (R)]I_LIEJ 
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(9) 



where P is the usual Kane matrix elementl 
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p = -i—(s\ Px \x) = h 

m 



2mo' 



(10) 



expressed in terms of the parameter Ep (Ref. [25]) and k ± - k x + 
ik y . We have also used that (S\p x \X) = (Sie y \Y) = (S\p z \Z). 
Equation (H)) is the 8x8 Kane HamiltonianU describing the s- 
type conduction and p-type valence bands around the F point 
in zincblend compounds. Note that the diagonal elements in 
Hamiltonian (g) correspond to the eigenenergies e/o of Eq. (Q): 
£io = £20 = ("conduction-band states," defined as the zero 
of energy), £30 = £40 = £50 = £60 = ~E g ("heavy" and "light" 
hole bands), and £70 = £so = —E g - A g ("split-off" hole band). 
Here, 



3^ {xl dV d 
4mlc 2 dy dx 



dV d 



dx dy 



W) 



(ID 



is the "atomic' 1 
Fig. 



1(a) 



SO parameter defining the split-off gap; see 
which schematically shows the conduction and va- 
lence bands of a zincblend structure. The circles indicate the 
k — eigenenergies. 

The Kane model treats exactly the conduction-valence band 
couplings within the truncated set of eight band-edge wave 
functions. It is important to emphasize that we have neglected 
contributions from the ^-dependent SOrterm in Eq. (Q), when 
constructing the Kane Hamiltonian (^).LI The SO interaction 
is accounted for only within the zone center Schrodinger's 
equation (pi) (parameter A g above). The diagonalization of 
the Kane Hamiltonian gives the dispersions E, h k around the T 
point. It is known that the Kane model presented here is not 
accurate for valence bandsU {e.g., wrong sign of the heavy 
hole masses). However, it provides a simplified and accurate 
description for the conduction electrons, which is the focus of 
our work. Next we discuss the Kane model in the context of 
heterostructures. 



2.3. Kane model for quantum wells 

Following Refs. [_] and |2(] we can straightforwardly gen- 
eralize the bulk Kane model of Sec. 2.2 to heterostructures. 



r 



Essentially, we have to introduce position-dependent (growth 
direction) band gaps which represent the diffe rent compounds 
comprising the heterostructure, e.g., Fig. Kb) . In this case, the 
form of the resulting Kane Hamiltonian is similar to that of 
bulk but with z-dependent diagonal matrix elements and with 
k z — > -id/dz. More specifically, defining E$ = H\i = H22, 
E& = H33 = H44 - H55 - H(&, and £7 = H77 = H^, we have 
for the double quantum well of Fig. 1(b) 

h 2 k 2 
E 6 = -— + V H (z) + h 6 (z), (12) 

Zmo 

E % = — + V H (z)-H(z)-E g , (13) 

2ntQ 

H 2 k 2 
Ei = — + V H (z) - h 7 (z) -E g - A g , (14) 

ZlflQ 

with k 2 - k 2 + k 2 , Eg and A g being the fundamental and split- 
off band gaps in the well region, respectively, and 



hi(z) = 6i h w (z) + 5bi h(z), 



i = 6, 7, ! 



(15) 



where h w (z) is a dimensionless profile function describi ng th e 



shape of a square well of width L w [and unit depth, Fig. |l(b)[ ] ; 
similarly , hb(z ) describes the shape of the central square bar- 
rier, Fig 



1(b) 



The parameters 6&, 67, 5&, Si,6, 5bi, and 6bs de- 
note the relevant band offsets between the well and the lateral 
and central barriers for conduction and valence bands. Defin- 
ing the z ero o f energy at the bottom of the conduction well 
[see Fig. 1(b)], we have 



6g = E w - Eg - 6 6 , 67 =6g+ A,, - A g , 
Sbg - Eb - E g - 6b6, Sbi - Sbs + Ai, - Ag. 



(16) 
(17) 



The corresponding expressions for a single well can readily 
be obtained from the above by setting the dbi's, to zero (i.e., no 
central barrier). 

Finally, note that we have added a "Hartree" potential Vh{z) 
in the diagonal elements; see Eqs. (|12|)-([l4|). The Hartree po- 
tential Vh(z) here contains contributions from the electron- 
electron interaction (mean field) relevant in quantum wells 
containing many electrons, the external gate potentials, and 



the modulation-doped potential (i.e., ionized impurities out- 
side the well region). In Appendix |a|, we describe in detail 
these distinct contributions to Vh and how they are calculated 
in our system. As we will see next, both Vh{z) and the struc- 
tural confining potentials contribute to the effective SO cou- 
pling for electrons. 
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(a) Energy band structure. 
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(b) Band oflfsets. 



Figure 1 : (a) Schematic of the band structure of direct gap zincblend 
semiconductors near the T point (k = 0). The label T(, represents 
the .? states in the conduction band, while T 7 (split-off holes) and T 8 
(heavy holes and light holes) represent the p-states in the valence 
bands, (b) Band offsets for a double quantum well of width L w with 
a central barrier of width L b . The relationships among the several 
offset parameters are given in (|16|)-(|l7[). 



3. EFFECTIVE SPIN-ORBIT HAMILTONIAN FOR 
ELECTRONS 

3.1. Folding down 

Since we are interested in SO effects for the conduction 
electrons, here we derive an effective Hamiltonian for them. 
To this end, let us rewrite our 8x8 Hamiltonian [Eq. (ph] in 
the block form, 



Hs* 



H c H, 



H n , Hy 



(18) 



where H c is a 2 x 2 diagonal matrix in the sector T(, (conduc- 
tion band) with identical diagonal elements £g [(|12|)] and H v 
is a 6 x 6 diagonal matrix in the sectors Tg and Tj (valence 
bands) with diagonal elements Eg [(|l3|)] and £7 [(J14|)], re- 
spectively. The 2x6 matrix H cv can be read off directly from 
the corresponding 2x6 block in Eq. (ph. 

Using the block form of our Hamiltonian (|l8|)) the eigen- 
value problem can be written in the compact form 



H c H„ 
H cv H c 






^v 



= E 



( \ 



(19) 



where ip c is a two-component spinor (conduction sector) and 
ip v is a six-component spinor (valence sector). Straightforward 
manipulations!— LI yield the effective Schrodinger's equation 



r H(E)ifi c = Eip c , 



(20) 



with 



<H(E) = H. + H CV (E - H.y'Hl 



(21) 



Table I: Truncated set of zone center wave functions u m (for simplic- 
ity we denote them by w/) used in constructing the matrix Hamilto- 
nian (|) 
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and £_ is a properly renormalized conduction-electron 
spinor£j 

The matrix elements of 'H(E) are given bytll 



•K(£)n = <H(Eh2 = E b + —{kl ri + L yi k : ), 
<H(E) n = W)\ x = ^k.[y 2 ,k z } = ~k-k z y 2 , 



where k\ — k ± k T —k 2 x + k\, and 

2 



E - Eg E - Ej 



-Ei\ 



n(z) = ( 

n(z) = ish; ~ e^e} 



(22) 
(23) 

(24) 
(25) 



We should emphasize that Eq. (|20|) is not really an eigenvalue 
equ ation as "HiE) depends on E. However, as we show in 
Sec. 3.2, we can still obtain a true eigenvalue problem by per- 
forming suitable expansions. 



3.2. Energy denominator expansions 



3.3. Projection into the quantum well subbands 



Since E g and E g + A are the largest energy scales in our 
system, i.e., 



X& 



Xi 



2 1,2 



2rriQ 



+ V H (z) - h(z) 



2 1,2 



h z k 

2niQ 



+ V H (z) - h 7 (z) 



E g + A g 



«1, 



« 1, 



(26) 



(27) 



we can expand the energy denominators in the y,'s [Eqs. ([2 
and (EST)] in the form 



y\ 



— (1-AT8 + --0 + 



l 



E g + A g 



1 n ^ l 

72 = — ( 1 — Yf, + ■ ■ ■ ) 

r E Q K A8 ' E„ + A t 



(1-AT7 + - ••) 
(1-AT7 + ---)- 



(28) 
(29) 



For the diagonal matrix elements < H(E)\\ - < H{E)22 we keep 
only zeroth-order {i.e., energy-independent) terms, while for 
the off-diagonal matrix elements < H{E)\2 = ( H(E) 2 , we keep 
in addition the first-order terms as they give the lowest non- 
vanishing contribution (because the off-diagonal matrix el- 
ements contain derivatives with respect to z). Straightfor- 
wardly, we then obtain the energy-independent one electron 
Hamiltonian 



<H(E) = H QW 1 + tj(z) 




ik + 



-ik. 





(30) 



where 



2 1,2 



Hqw - 



h z k 



n 2 a 2 



2m* 2m* dz 2 



V sc (z), 



(31) 



[the subscript "sc" emphasizes that the potential is to be.de- 
termined self-consistently - see Appendices |a| and|B[] andllJ 

3n 2 \E ! , Ee + AJ- 



1 1 2P 2 I 2 



m* m Q J(1 ^ g ^g^^g, 
V sc (z) = V H (z) + S 6 h w (z) + 6 bb h b {z), 

dh w (z) dh b {z) dV H (z) 

T](Z) = T] w -J^- + Tlb—T_ r lH- 



dz 



dz 



dz 



with 



p 2 

** = T 

p 2 



i 
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L£f ~ (Eg + A g ) 2 

E\ (Eg + Ag) 2 

6bx _ 5bi 
El (Eg + Ag) 2 



(32) 
(33) 
(34) 

(35) 
(36) 
(37) 



Here we define a quasi-two-dimensional (2D) model start- 
ing from the three-dimensional (3D) Hamiltonian (p0|). The 
idea is essentially to obtain a 2D effective model similar to the 
well-known Rashba model, but now for the case of wells with 
many subbands. To this end we (i) first project [(BOb into the 
spin-degenerate eigenstates of Hqw [Eo_ (plj)] (note that Hqw 
does not contain the SO interaction) :Ej Iknv)^ = |k M v)|cr z ), 
<r|k„v) = exp(/k ll -r ll )9?,,(z), v = 0, 1,2, . . , and cr z = + (or f , J,), 

which correspond to the subband energies £j l|V = -^ + S v , 
with fi v being the quantized levels of the well, and then (ii) 
consider a reduced set of subbands (e.g., two) by truncating 
the basis set used. In this section we simply assume that we 
know the eigensolutions of Hqw', later on we actually calcu- 
late them within a self-consistent procedure, from which we 
can explicitly determine the relevant SO coupling constants in 
our problem. 

The matrix elements ofJ-{(s) in the {|k M v) (7 -J basis are 



<k„vK±|*WCE)|k„v' 
(k ll v|(±[-K(£)|k ll v')|+) = +/^^, 



x ,h 2 kl 

2m 



with the generalized SO couplings 

n , - T H + r"' + T b 



where 



-Vh(v\ — ^-—\v), 



r*. = +vb(v\ 



dz 
dh w (z). , 



dz 
dh b (z) , 



Iv'>, 



dz 



Iv')- 



(38) 
(39) 

(40) 

(41) 
(42) 
(43) 



The coefficients T^,,, F*,, and T*,, denote the contributions 
from the Hartree potential, the quantum-well edges, and the 
central barrier edges, respectively. It is convenient to split the 
Hartree contribution into two terms, i.e., Vh(z) = V e (z)+V g (z), 
where V e (z) is the purely electronic Hartree potential and 
V g (z) denotes the contributions from the external gate po- 
tential and the modulation doping potential. Hence T V J - 

-t]h(v\ 4. \v') - t]h{v\ j. \v'). This separation will be use- 
ful when discussing our results. 

We emphasize that the diagonal (in v, v') parameters r\ vv 
correspond to the Rashba coupling in the vth subband, i.e., 
a v - 7w. The off-diagonal terms 7] vv > arise due to the inter- 
subband coupling. Interestingly, these new SO terms can be 
non-zero even in structurally symmetric wells, since they arise 
from quantum-well states of distinctive parities. 

For completeness we present here the linearized Dressel- 
haus couplingsLJ in the vth subband 



P v =Mv\kj\v), 



(44) 



where the constant /3d is the bulk Dresselhaus SO parameter. 



D 



We can easily rewrite the above expression in the more con- 
venient form 



2/72* 

Pv=Pd-^[&v-(v\V{z)\v)], 



(45) 



In Sec. mwe shall use the above form to discuss how the Dres- 
selhaus couplings vary as a function of the system parameters. 



where 



e ± (n) = ^W«) 2 + (6 t|l _ ± a-fc,,) 2 , e ±ie -- 



h. 



3.4.3. SO-induced effective mass renormalization 



3.4. Two-subband case 



To illustrate the procedure of Sec. 3.3 , let us explicitly work 
out here the case of a quantum well with only two subbands 
v = 0, 1 . In Sec. |5| we shall investigate in detail the SO 
couplings for single and double quantum wells with two sub- 
bands. 



3.4.1. 4x4 Hamiltonian 

With the basis ordering {|k„0) T , |k„0)j., |k„l) T , |k,,l>j,}, 
Eqs. (f38|) and (J3^) yield the effective Hamiltonian 



£*„o -iaok- -vqk- 



H = 



ia k+ & ku o ii]k + 

-irjk- £fc n i -ia\k 
irjk+ iat\ k+ S^i , 



(46) 



where the Rashba couplings are giijen by a v = t] vv , v = 0, 1, 
and the intersubband SO couplingUby r\ = 7701 [see Eqs. ( |40| ) 
and @] and 



G k „v = o v + v = 0, 1. 

2m 



5.4.2. Eigensolutions 



(47) 



The energy eigenvalues S^ of Eq. (|-6|) are straightforward 
to obtain: 



£h,M,M = £ *,+ + A 2 a + k « + A i yji^ku) 2 + (£*„- + A 2 a-k ll ) 2 , 

(48) 
where /I2 = + are spin quantum numbers and A\ — + are the 
subband (or pseudo spin) indices, and 



£>h± = »(£*,, 1 ±£jt„o). oc± - -(ai + ao). 
The corresponding (normalized) eigenvectors are 

1 -iAiT]k u e- w ^ 



(49) 



Hi,A 2 >= Jl+^ 



Q 2 (0) 



^(1/) 



A\A.2t]k n 



e Al {ji) + Me M (0) 
-iA 2 e- ie 

1 



»iku-r,. 



4/T 



(50) 



Expanding the energy dispersions |Eq. (| |)| around A M = 0. 
we obtain to second order 

h 2 k 2 
£>k u ^o,A u ^ ~ <5+ + A\ &- + A 2 (a+ + A { a-)k„ + — -£-, (52) 

Zm ' 

M 

where m\ are the effective masses 



(53) 



1 ± 



2£.so 
~AG 



where S so = \m*T] 2 /h 2 and AS = 2£L. Note that the mass 
renormalization is solely due to the intersubband-induced 
SO coupling 77. For the realistic wells we investigate here 
2S S0 /AS « 1 for single wells but can reach ~ 0.1 for double 
wells (Sees. I and |). 



3.4.4. Determining the SO couplings 

As mentioned previously, we determine the SO orbit cou- 
plings (here specifically ao, a\ , and 77) from the self-consistent 
eigensolutions of the quantum well without spin orbitJJ via 
Eqs. (HQ))-(p3[))- In Sec. [4] we detail the quantum well sys- 
tems investigated and briefly outline the self-consistent pro- 
cedure used to obtain the eigensolutions (a full description is 
provided in Appendices M and B). We then present results for 
single and double wells with two subbands; i.e., we calculate 
ffo, a\, and 77 and discuss in detail the several distinct contri- 
butions to each of these quantities. 



4. QUANTUM- WELL SYSTEMS AND 
SELF-CONSISTENCY 

Figure g shows a schematic view of the quantum-well sys- 
tem we study: a well of width L„, centered at z — (z: growth 
direction), and two adjacent symmetrically doped regions of 
widths w in the barriers. We also consider double wells by 
inserting an additional (central) barrier of width Lj, in the well 
region. The doping densities of the left and right regions, p„ 
and pi,, respectively, can be used to control the degree of struc- 
tural inversion asymmetry of the wells (in Sec. pi however, we 
present results only for p a - pi,)- The external gates V a and 
Vb, located at the end points ±L, can also be used to control 
the degree of inversion asymmetry and to vary the areal elec- 
tron density in the well. 

Since our wells have many electrons and are subject to ex- 
ternal gates, we have to solve the Schrodinger and Poisson 



equations self-consistently ("Hartree approximation"!— I) in or- 
der to determine their potential profile V sc (z) [see Eq. Q33J>] 
and corresponding eigenfunctions and eigenenergies. In Ap- 
pendices |a| and |b] we describe in detail our standard self- 
consistent procedure. 
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Figure 2: Schematic view of our quantum-well system. The doping 
densities p a , p b and the external gate voltages V a and V b can be used 
to control the degree of the structural inversion asymmetry. 

Before going into the discussion of the SO couplings in 
detail, let us first have a look at the outcome of a typical 
self-consistent simulation we perform. Figure 3(a) shows 
the self-consistent potential V sc (thick solid line) for a sin- 
gle well with two subbands; the corresponding self-consistent 
wave functions 4>o(z) and i//i(z) are also shown. The energies 
of the two lowest subband edges (see levels in the well) are 
£o = 309.09 meV and & x = 406.39 meV (A«5 = 97.3 meV). 
Here we fix the chemical potential at p — 413.40 meV with 
respect to the V — origin ("constant chemical potential 
model", see below) and set the external gates to V a = and 
Vb = 1200 meV. The two subbands are occupied with areal 
densities « = 18.7422xlO n cm 2 and^ = 1.2578x10" cm 2 , 
respectively. The electronic Hartree potential V e (short dashed 
line) and the the external gates (plus modulation doping) con- 
tribution V g (long dashed line) are also shown. Figure 3(b) 
shows the corresponding "force fields" F e - -dV e /dz arising 
from the confined electrons in the well and F g = -dV g /dz 
coming from the doping regions (±12 nm to ±18 nm) plus 
the external gates (F g and F e will be useful when discussing 
the SO couplings further below). Using the self-consistent 

tforwardly 

©-(Hi] 

-3.81 meVnm, oq - -5.44 meVnm, a\ - -3.74 meVnm 
ySo = 0.87 meV nm, and/3i = 2.50 meV nm. 



solutions 4i v {z), v - 0, 1, we can straightforwardly calcu- 
late the relevant SO couplings [via Eqs. (Q)-(^3|)]: r\ - 
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Figure 3: (a): Self-consistent potential energy V sc (thick solid line) 
and the corresponding wave functions tf/o and tf/[ for the single well 
Alo.48lno.51As/Gao 47 Ino.53 As with external gates V a = Oe Vb = l.2eV 
(see Fig. Q). The the electronic Hartree potential V e (short dashed 
line), the external gate plus modulation doping contributions V g (long 
dashed line), and the corresponding force fields F e = -dV e /dz and 
F g = —dVg/dz are also shown in (b). The two levels in the well (solid 
lines) denote the energies of the first and second subband edges, 
while the dotted level indicates the chemical potential. 



5.1. Single wells 



5. RESULTS 



5.1.1. Single-well parameters 



Here we focus on single and double quantum wells with 
only two subbands. More specifically, we calculate three SO 
couplings: the intersubband-induced SO coupling 77 = 7701 and 
the two Rashba-type couplings ao = Voo an d a \ = Vu- We 
consider two experimentally relevant cases: the constant areal 
density («7--constant) and the constant chemical potential (p- 
constant) models. In our simulations we always keep V a = 
as a reference potential and vary V&; see Fig. For complete- 
ness, we also calculate the two Dresselhaus constants /3o and 
fi\ [see Eq. (44)] within each subband. 



We conside 
quantum well 



realistic Alo.4sIno.52As/Gao.47Ino.53As single 
We assume doping densities p a - Pb — 



4 x 10 18 crrT 3 with widths w — 6 nm ("sample 3" in Re£ |33p . 
Table |Ey summarizes band parameters, potential offsetsJJ well 
widths, and other important parameters of our system. The 
coefficients 77,,. and r/ H in Tab. y are defined in Eq. ( pq ) and 
), respectively. Here, the Dresselhaus parameter /3D in 
) is assumed to be the same as that of the GaAs (see 




Table II: Relevant parametersLj (see Fig. |lj) (at 0.3 K) for the sin- 
gle quantum well Alo.48 Ino.52As/Gao.47lno.53As system in our simu- 
lations. The doping regions have widths w = 6 nm and densities 
Pa = Pb = 4 x 10 18 cnT 3 (see Fig. 0). All energies are in eV and 
lengths in nm; the coefficient p™ is in nm 2 and ij w in meV nm 2 . The 
Dresselhaus coupling constantU (3D is in meV nm 3 . 
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5.1.2. SO couplings: single wells 

Figure [| shows the strength of the Rashba (a v , v - 0, 1 ; 
dashed lines), Dresselhaus (J3 V , v = 0,1; dotted lines) and 
intersubband-induced (77, solid line) SO couplings as func- 
tions of the gate voltage Vi, for bot h the nj-constant and 
the //-constant models, Figs. |4(a)| and 1(b), respectively. At 
Vb - V fl = eV, our sample is completely symmetric and, 
as expected, the Rashba couplings ao and a\ are zero. We 
note that the Dresselhaus couplings /Jo and fi\ are practically 
constant in both models. This follows from Eq. (W5t) which 
shows that in each subband the Dresselhaus coupling is essen- 
tially the difference between the expected value of the self- 
consistent potential in the respective subband and the corre- 
sponding eigenenergy. The Rashba couplings, on the other 
hand, vary considerably with Vb, although showing a similar 
trend in both models. Interestingly, they change signs about 
Vb - (symmetric configuration), but always with |ap | > |o fi|. 
Our calculated ao within the yU-constant model [Fig. 4(b) ] is 
consistent with the measurements of this quantity by Koga et 
nl I I I whose samples have a constant chemical potential. 

The new i ntersu bban d-indu ced coupling 77 [see the solid 
lines in Figs. 4(a) and 4(b)] is non-zero even in the sym- 
metric well configuration (Vj = = V a ). It has a strength 
comparable to the Rashba and is at least twice as large as the 
Dresselhaus. In contrast to the Rashba couplings, the intersub- 
band SO 77 does not change sign with Vb- In fact, for the single 
well investigated here 77 is almost constant with Vb, although 
it varies slightly mo re in the //-constant model [compare the 
solid curves in Figs. 4(a) and 4(b)]. 

To more easily understand the results above, we analyze the 
several contributions to the SO couplings separately. To this 
end, we rewrite [see comments following Eq. (S3)] 77,,,/ for a 
single well in the form 



77T 

'/IT' 



= r e , + r* 

A vv' v 



+ r M 



(54) 



where we have set F'', = in (40), i.e., no central barrier 



contribution, and have split the Hartree contribution into its 
purely electronic HL and the external gate (plus doping po- 
tential) I*, parts. Hence, for two subbands, each of the SO 
couplings has three contributions: 77 = 77^ = T e m + T g Qy + VZ V 
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Figure 4: Rashba a, Dresselhaus /J and intersubband-induced rj SO 
coupling constants for the Alo.4gIrio.52As/Gao.4fIno.53As quantum well 
as functions of the gate voltage Vb (see Fig. g). In (a) the total 2D 
electron density is kept constant at n T = 20 x 10" cirT 2 and in (b) 
the chemical potential is kept constant at p = 200 meV. 
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Figures 5(a) - 5(c) show the above contributions separately for 
the fty-constant case (similar results hold for the yU-constant 
model, in the parameter range studied). 



Figure 5(a) shows that the external gates and doping con- 
tributions to 77 (r£j curve) are essentially zero, while the elec- 
tronic Hartree contribution (T1, curve) and the structural (F^ 
curve) contributions are comparable in magnitudes and both 
negative. In contrast, for both ao and a\ the largest con- 
tributions come from the external ga tes to gether with dop- 
ing reg ions [see the curve Fq in Fig. 5(b) and the curve T s n 
in Fig. 5(c)]; these account for 60% of ao and 100% of a\. 
The electro nic H artree contribution is negligible in ao [curve 
Tqq in Fig. 5(c)] while the structural part (TJJ,) accounts for 
about 30% of it. On the other hand, the structural and elec- 
tronic Hartree contributions in a\ essentially cancel out (same 
magnitude and opposite signs); cf. the T^ and T^ curves in 
Fig-|5(4 
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Figure 5: Several distinct contributions to the coupling constants 77 
(a), a (b) a nd ai (c) for the single GalnAs quantum well shown 
in Fig. 4(a) (n^-constant model) as functions of the external gate 
V* (Va = 0). These contributions arise from: the electron density 
(Hartree potential), the external gate (together with donor regions), 
and the structural well potential; these are denoted by the superscripts 
e, g and w, respectively. 



Rashba couplings a v are strongly m odul ated by external gates. 
This is even more so for a\, Fig. 5(c) , for which the struc- 

Looking at the 
-dV e ldz in 



tural and electronic contributions cancel out. 

wav e fun ctions \J/q and \jj\ and the force field F e 

Fig. 3(b) , we can see that the electronic Hartree contribution 

( \v\Fe\v)) is almost zero (though slightly negative) for the 

lowest subband and positive for the first subband. The struc- 
tural well contributions F"j, [see Eq. (|42|)] to a v are similar for 
both subbands, though |njJJ > \T K n \, because the nonzero bi- 
ases (Vi, + 0) cause the wave functions to shift toward one 



side of the well [e.g., Vb — 1.2 eV in Fig. 3(a)] 



On the other hand, the contribution tT. ■ -(0|F g |l) to 
the intersubband coupling 77 i s ess entially zero since the wave 
functions [i/'o and \j/\ in Fig. 3(a) ] are orthogonal and, again, 
F g is constant. Hence 77 is not as sensitive to the external 
gates as the Rashba couplings. Most of the modulation of 
rj arises from the electronic Hartree and structural contribu- 
tions, which both have the same sign and magnitude as shown 
in Fig. 5(a). 



5.2. Double well 



5.2.1. Double-well parameters 



Table 



111 



shows the band parametersciltJ for the dou- 
ble quantum well Alo.4Ino.6Sb/InSb with one central barrier 
InSb/Alo.nlno.gsSb. Hereafter we refer to this heterostructure 
as InSb double well. The meaning of some of these parame- 
ters (e.g., band offsets) can be seen in Figs, p] and 0. 



Table III: Relevant parametersUU (at 1 K) for the InSb double well 
see Figs, m and Eh. The width of the doping regions is w = 4 nm and 
their densities are p„ = p b = 3 x 10 18 cm 4 . All energies are in eV 
and lengths in nm. The coefficient n H is measured in nm 2 while 77,,. 
and rj b are measured in meVnm 2 . The Dresselhaus constant (3D is 
measured in meV nm 3 . The parameters in the last column are to the 
InSb binary compound. 
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in Fig. 3(b) 



We can understand the above remarks by looking at the self- 
consistent potentials and the "force fields" F e - -dV e (z)/dz 
(short dashed curve) and F g = -dV g (z)/dz (long d ashed 
curve) - note that T' , ~ (v|F,|v) 
This figure was obtained for Vb 
play the general behavior for all quantities shown. The force 
field F g is essentially constant, except within the donor re- 
gions where the wave functions are vanishinglysmall. Hence, 
the matrix element (v|F g |v) [see Eqs. Q4C|) — (H-3[)] is approxi- 
mately linear in the external gate Vb- This explains why the 



5.2.2. SO couplings: double-well case 

Figure ^ shows the Rashba a, Dresselhaus /3, and 
intersubband-induced 77 SO couplings as functions of the gate 
voltage Vi, (here again V a = 0) for both the n^-constant (a) 
and jU-const ant ( b) models. We first discuss the 777-constant 
model [Fig. 6(a) ], Here the Rashba couplings (dashed lines) 
are most sensitive to the external bias Vb, being essentially 



10 



the largest of all SO couplings for very asymmetric structures 
(i.e., high biases). The Dresselhaus couplings (dotted lines) 
are almost identical (J3q ~ /3\) and mostly independent of the 



Figs. 7(a) and 7(b)]. 
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Figure 6: Rashba a, Dresselhaus /3 and the intersubband-induced t] 
SO couplings for a InSb double quantum well as functions of the 
right gate voltage Vb- In (a) the total electron density is kept constant 
at n T = 10 x 10" cm -2 and in (b) the chemical potential is kept 
constant at// = 100 meV (relative to to initial bottom well). 



external gates. The SO coupling r\ (solid line) is an even func- 
tion of the external gate Vi, and presents a "resonant behavior" 
around the Vb — V a - eV configuration, at which our sample 
is symmetric. While the Rashba couplings are both zero at this 
symmetric configuration, we note that they are odd functions 
of the external gate (with |aol > [ail), have opposite signs and 
abruptly change magnitudes around Vj, = ( over a 40 meV 
wide region). For the /j-constant model [Fig. |5(b)| ], a similar 
picture as above also holds; note, however, that in contrast to 
the Mr-constant model, in the /i-constant case the positive and 
negative bias configurations are not equivalent as they corre- 
spond to the well having different numbers of electrons. 

For completeness we show in Fig. M the behavior of all cou- 
pling constants near the sym metric point Vb — V a — eV for 
the double well in Fig. 6(a). Note that the Dresselhaus cou- 



plings /Jo andySi pres ent a (double) crossing over a 160 meV 
wide region [see Fig. 7(b)]. However, this is a minor effect: 
note the change in the scale of the vertical axis. While the 
resonant behavior of 77 is accompa nied by an enhancement of 
about 10 in its magnitude [see Fig. 5(a)], we see no substantial 
change in the magnitudes of the/3's near the zero-bias case [cf. 
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Figure 7: Rashba a, intersubband-induced rj and Dresselhaus /3 cou- 
plings vs Vb about the sym metric configuration Vi, = V u = eV for 



the double InSb well in Fig. 6(a) 



The relative strengths of the Rashba and Dresselhaus cou- 
pling constants to the intersubband-induced SO coupling are 
shown in Fig. pi The Rashba couplings have the largest 
strengths (note the pre-factors in front of a v /rj in the legends). 
In contrast to j3 v /t], the linear behavior of the Rashba ratios 
a v jr\ near Vb - (see insets) shows that a v and 77 undergo 
similar variations near the symmetric configuration. As ob- 
served before, the intersubband-induced coupling r\ becomes 
important near Vi, — (Fig. [|). 
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Figure 8: Ratios a v jr] and f3 v lr] for the InSb double well in Fig. W. 
The insets are blowups around Vb = 0, the symmetric configuration. 

Figure g (similar to Fig. g for the single-well case) shows 
the several contributions to each of the SO couplings 77, co 
and a\ for the double-well case. Here, in addition to the elec- 
tronic Hartree, the gate (+ doping regions), and the well con- 
tributions, there is an additional structural term arising from 
the centra l barrier (superscript b). A general feature in Figs. 
9(a)| -9(c) is that the structural contributions (well and central 
barrier) almost cancel out because they have opposite signs 
(see the curves with superscripts w and b). These terms have 
opposite signs because the derivatives dh w (z)/dz (well) and 
dhb(z)/dz (barrier), which enter the coupling constants [see 
Eqs. (|40|), (|42|), and Q], have opposite slopes. Similarly to 
the single-well case [Fig. 5(a)] the contribution of the exter- 
nal gates (which includes the doping regions) to the intersub- 
ban d SO coupling 77 is vanishingly small [see the F£, curve in 



Fig. 9(a)]. Hence, 77 is mostl y due to the electronic Hartree 
contribution [curve F e m in Fig. 9(a) ]. In addition, the gate con- 



tribution to ffo and a\ for the InSb double well is linear in Vb 
as for the single-well case. Hence, the Rashba couplings a^ 
and a\ for the double InSb well are essentially determined by 
the electronic (Hartree) contribution and are modulated by the 
gate contribution. Summarizing: looking at Fig. pi we can see 
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by looking more closely at the self-consistent wave functions 
of the InSb double well around the symmetric configuration 
(Vh - 0). The top row in Fig. fTO shows the self-consistent 
potential profile of the double well and the normalized wave 
functions i/>o (short dashed line) and \j/\ (long dashed line) for 
the lowest v = and for the first excited v = 1 subbands at 
three distinct gate voltages: V/, = +0.3 eV, Vj, = eV and 
Vb - -0.3 eV (left, center, and right columns, respectively). 
For positive bias i/>o is mostly localized in the left well and i//\ 
in the right well, while for negative biases this configuration is 
reversed. The electronic Hartree contribution to the potential 
energy V e and the corresponding force field F e = -dV e /dz are 
shown on the second row, thin and thick lines respectively. 
Notice that F e is practically zero in the central barrier re- 
gion (-10 < z < 10 nm) and has opposite signs within the 
wells (-25 < z < -10 nm and 10 < z < 25 nm). Hence 
the quantities F°°(z) = ifr (z)F e Mz), F] l (z) = Mz)F e >fn(z) 
and F® l (z) = iffo(z)F e ij/\(z) have the forms shown on the third 
and fourth rows. The integral over z of these quantities de- 
fines the electronic Hartree contributions to the spin-orbit cou- 
plings aso, a\, and 77, i.e., r e m ~ (0\F e \Q), F n ~ (l[F e |l), and 
(0|/v|l), respectively. Since the electronicHartree con- 



tributions dominate over the others, see Figs. |9(a) -9(c), the 
abrupt changes in the Rashba couplings and the resonant be- 
havior of 77 around Vb - follow straightforwardly. 
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Figure 9: Different contributions to the SO couplings 77, ao and ct{ 
for t he I nSb double well in the constant areal density model [see 
Fig. 5(a)] as functions of the external gate Vt (V„ = 0). In the sub- 



figures we show the contributions to the coupling constants coming 
from the areal electronic density, indicated by the superscript e, and 
from the external gate + donor regions g, and also from the structural 
potential, being w for the well and b for the central barrier. 



that (i) the structural contributions (well and barrier; dashed 
curves) almost cancel out, and (ii) the external gate (dotted 
curves) modulates the Rashba couplings or v ; therefore, for the 
double well investigated here (iii) most of the strength of these 
three coupling constants (77, ao, and a\) comes from the elec- 
tronic contribution (dot-dashed curves). 

It is instructive to investigate in more detail how the reso- 
nant behavior in r\ comes about, as well as the abrupt changes 
in the Rashba couplings; see Fig. H. This can be accomplished 



5.2.3. Density anticrossings and effective masses 

Figures |lT](a) and [H|(b) show anti-crossings of the areal 
densities nj for thpJnSb double well near the symmetric con- 
figuration Vb - 0l_i where the strength of the intersubband- 
induced SO coupling 77 is the strongest (-16.7482 meVnm) 
while the the energy difference between the subband edges 
AS - Si - So (0.9353 meV) is the smallest. In accord 
with Eq. (p2|), we find an appreciable change in the bulk 
effective mass m* near k„ -■- Olj. The ratio S W /AS [see 



Eq. (|53h] is shown in Fig. |1 l|(c) and the ratio m ± - m* ± /m* 
in FigTTTTVd). These intersubband-SO-induced changes in the 
effective masses m* ± may have a sizable effect on the measured 
mobilities and cyclotron frequencies in InSb wells. 



6. SUMMARY 

Starting from the 8x8 Kane model in heterostructures, we 
have derived in some detail an effective electron Hamiltonian 
which contains a new intersubband-induced SO interaction 
term which arises in quantum wells with more than one quan- 
tized subband. Unlike the usual Rashba SO term, the intersub- 
band SO coupling here is non-zero even for symmetric wells. 
For structurally asymmetric wells we have also accounted for 
the Rashba-type SO interaction within each subband. 

We have also outlined the projection procedure ("folding 
down") to obtain quasi-2D Hamiltonians by integrating out 
the confined variables. For two subbands in asymmetric wells 
we find a 4 x 4 quasi-2D Hamiltonian resembling Rashba's, 
but containing three SO couplings: the two Rashba couplings 
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Figure 10: To p row : Wave functions i//q and t/o for the InSb dou- 
ble well (Sec. 5.2. lb with external biases V b = +0.3 eV, V b = eV 



and V b = -0.3 eV (1st, 2nd and 3rd columns, respectively). The 
electronic Hartree contribution V e (z) to the potential energy and the 
corresponding force field F e = -dV e /dz are shown on the second 
row. The third row shows the quantities Ff(z) = tffoF e if/ and 
F\ l (z) = i//iF e ij/i and the bottom row shows Fj? (z) = if/oF e tf/{. Here 
n T = 10 x 10 11 cm -2 . The vertical dotted lines marks the center of 
the 20 nm wide AllnSb barrier within the 50 nm wide InSb well. 



ffo and a\ and the intersubband SO coupling r\. For this two- 
subband case, we have investigated thoroughly these three SO 
couplings for realistic modulation-doped single and double 
wells. By performing a detailed self-consistent calculation in 
which we solve both Poisson's and Schrodinger's equation it- 
eratively, we have determined the strengths of ao, a\, and r\. 

Each of these coupling strengths contains contributions 
arising from (i) the potential-well (and barrier) offsets, (ii) the 
electronic Hartree potential, and (iii) the external gate poten- 
tial plus the modulation doping potential. We have performed 
our simulations by either keeping the areal electron density 
iir in the well fixed (/^-constant model) or by keeping the 
chemical potential ju fixed (/i-constant model). In the parame- 
ter range investigated, both models give similar result for the 
calculated SO couplings. 

For the single well investigated, ao is mostly determined 
by the external gate (+ modulation doping) contribution; with 
the structural + electronic Hartree being about half of that of 
the gate. On the other hand, a\ is essentially determined by 
the external gate (+ modulation doping) contribution, since 
the electronic Hartree and the structural contributions cancel 
out. The intersubband SO coupling r\ is essentially determined 
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Figure 11: (a) Total electronic areal density n T = no + ri\, (b) sub- 
band areal densities n andni, (c) ratio £ so /A£ and (d) effective mass 
ratios m± = m* ± /m* [Eq. (p3J)], as functions of the external gate V b 
(V a = 0) in the //-constant model (ji = 100 meV). 



by the electronic Hartree potential and the structural potential 
contributions (both of the same order); the external gate (+ 
modulation doping) potential contribution to rj is nearly zero. 
Hence, while ao and a\ can be modulated by the external gate 
potential, r\ is only slightly influenced by it. 

For double wells the SO couplings show more peculiar 
behaviors. While the Rashba couplings »o and a\ abruptly 
change magnitudes and signs around the symmetric configu- 
ration (zero external bias), the intersubband-induced SO cou- 
pling presents a resonant behavior being enhanced by a factor 
of 10 (with no sign change) around this point. For the double 
well investigated the structural contributions to ao, a\ and 77, 
due to the potential offsets of the edges of the well and the 
central barrier, cancel out. In addition, the contribution of the 
external gate (+ doping region) to 7/ is vanishingly small (as 
for the single-well case). Interestingly, the dominant contri- 
bution to all three SO couplings ao, a\, and r\ comes from 
the electronic Hartree potential. However, this contribution is 
highly influenced by the external gate, particularly around the 
symmetric configuration as the electrons can easily localize in 
either well for slight (positive or negative) changes in the gate. 

Finally, we have also calculated the effective mass renor- 
malization due to the intersubband SO interaction (the 
Rashba-type interaction does not produce a mass change). For 
the double well investigated, we find that this mass renormal- 
ization is the largest (~ 10%) around the symmetric poten- 
tial configuration (zero external bias), for which the splitting 
of the two subbands is the smallest. This mass change can 
possibly have an effect on mobility and cyclotron-resonance 
measurements. 
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Appendix A: SELF-CONSISTENT PROCEDURE 



b. Numerics 



1. Effective Schrodinger equation 

The single-particle electron Hamiltonian Hq W of our quan- 
tum wells [Eq. (|3l|))] is clearly separable. The transverse mo- 
tion {x,y) is free while that along the z direction is confined 
by the quantum well. To solve the corresponding Schrodinger 
equation HQ^P^fr) = E^W^^r) we assume a wave func- 
tion of the form 



¥ k|iV (r) = (r|k„v> = -L exp(/k M ■ r n )^ v (z) (Al) 

VA 



(A is a normalizing area) which leads to the ID Schrodinger 
equation 



We use the sixth- 
Schrodinger equation.L 



merov method to solve the 
Poisson equation (see Sec. ||) 



is solved via a semi-analytical Numerov method.li!] All nu 
merical integrations are performed using a Gaussian integra- 
tion methodlij In our numerical implementation we use the 
dimensionless form of Eq. (A2), 






0V = ^v(z). 2 = 



/' 



where 



V v 



—[v sc (z)-e v ], 



nh 2 



e\ 



*P ' 



(A3) 



(A4) 



We choose / = 1 nm as our length unit and s\ as the relevant 
energy scale. 



■^ + V ° AZ) 



h 2 k 2 

2m* 



£ v ^ v (z), (A2) 



Poisson equations for the electronic and gate plus 
modulation doping potentials 



from which we obtain the quantized energy levels £„ and wave 
functions if/ v (z). As we shall see, the subband structure of the 
well Ek ]]V - & v + ft 2 k 2 /2m* and the corresponding total wave 
function H^ ,,(r) will be used (within a self-consistent proce- 
dure) to construct the electron charge density, from which the 
corresponding Hartree potential can be obtained via the Pois- 
son equation. 



As mentioned in Sec. 3 2 V sc in Eq. ( A2) contains not only 



the structural confining potential but also the "Hartree contri- 
butions" (i) the purely electronic mean-field potential (elec- 
tronic Hartree potential) and (ii) the external gate potential 
plus the modulation doping potential. Further down we dis- 
cuss these contributions in detail. Each of these contributions 
is determined from a Poisson equation with an appropriate 
charge distribution and boundary condition. 



a. Self-consistency 

Since the electronic charge distribution p e (z) 
(oc 2 V K r |k M v)| 2 ) depends on the detailed form of the sev- 
eral potentials (modulation doping, gates, and electronic 
Hartree), and these, in turn, depend on p e (z), we have to 
solve the problem self-consistently. The standard procedure 
is as follows: (i) to solve Eq. (A2)) with an initial guess for 
V sc which we take to be just the structural potential plus the 
external gates and modulation doping potential [i.e., in the 
first run we do not include the electronic Hartree potential 
Ve(z)]> (ii) to construct the electronic charge density p e (z) 
[from the eigenfunctions obtained in step (i)] and the corre- 
sponding V e (z) via Poisson equation; and (iii) to solve again 
the Schrodinger equation with the new V sc , which in this new 
iteration includes V e (z) (as well as the other potentials: gates, 
modulation doping, and structural confinement). We repeat 
this process until convergence is attained. 



The self-consistent electronic potential energy V sc (z) - 
-e4> sc (z) can be split in two parts, V sc (z) - V w b(z) + V H (z). 
V w b(z) - V w (z) + Vb(z) described the structural quantum-well 
potential. The "Hartree" contribution Vh(z) - V e (z) + V g (z) 
arises from the electronic charge density and from the exter- 
nal gates plus the modulation doping regions (symmetrically 
located around the well; see Fig. |T2|). Figure |l2| also shows 
the Dirichlet boundary conditions V a and V/,, which are in fact 
the external gates at the end points +L of our system. 
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Figure 12: Schematic representation of doping layers of width w and 
densities p a andp b plus the external gates V„ and V b . By varying the 
external gates (we usually keep V a = and vary V/,) we can alter the 
spatial symmetry of our quantum wells, Fig. B. The curves illustrate 
the calculated gate + modulation doping potential V g (z) for V a = V/, 
(dashed line) and V a > V b (solid line), both withp fl > p b . 



a. Gate+modulation doping potential 

We can write separate Poisson equations for V g (z) and V e (z) 
as these arise from distinct charge densities. For V g (z) we have 
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(see Fig. ph 



P> 



e r eo 



Pa, 

0, 

Pb, 

0. 



fO, (I) : -L < z < -L d , 

(II) : -L d < z < -L d + w, 
(III) : -L d + w<z<L d -w, (A5) 
(IV) : L d - w < z < L d , 
(W):L d <z< L, 

where eo is the permittivity, e r is the dielectric constantLrU 
andp„j, are the doping densities. From the continuity of V g {z) 
and its first derivative and assuming the Dirichlet boundary 
conditions V g (-L) = V a and V g (+L) = Vb, we find 

'c\z + c 2 , (I) : -L < z < -L d , 
\Az 2 + c 3 z + c 4 , (II) : -L d <z< -L d + w, 
V g - I c 5 z + C6, (HI) : -Ld + w <z< L d - w, (A6) 

\Bz 2 + c 7 z + c 8 , (IV) : L d - w < z < L d , 
c 9 z + c m , (V) :L d <z<L, 



with 



e 2 p a 



Pb 



(A7) 



where the constants c,- are given in Appendix 0. Figure [T3 
shows two solutions of Eq. (A5), both having p a > pi, and 
V a - Vb (dashed line) and V a > Vb. 



b. Electronic Hartree potential 

The electronic Hartree contribution V e (z) is determined 
from 



j2 V e<Z) = ~Pe(z), 



e 
e r eo 



with (including spin) 

2e 
A 



2e \-i t em* 

Pe(z) = — Jj \^(z)\ 2 f(E k J = -^ k B TA e (z), 

hi 



(A8) 



(A9) 



where 



Uz) = J] l-A„(z)l 2 ln[l + e/XrtVW], (A10) 



and 



/(£*„v) = 



l 



2 7,2 



1 + s (E hll ,-fi)/k B T ' 



Ek„v = 



h z k; 
2m* 



+ G V . (All) 



We sobie Eq. (A8) for V e (z) using an accurate Numerov 
schemeU with the Dirichlet boundary conditio ns V e (+L) - 0. 
Similarly to the Schrodinger equation in Eq. ( A3 )), we find 
it convenient here to write the Poisson equation (A8) in a di- 
mensionless form 



d 1 - 



/I, 



ksT 



lAAX), s 2 V e = V e (z), (A12) 



where ej is the energy scale given in Eq. ( A4 ) and 

e 

£2 



e r eol 



(A13) 



c. Electron density and chemical potential 

From the total electronic charge 

dVp e (z) = en T A (A14) 



/• 



we can straightforwardly [using Eq. (|A9j))] obtain the total 
areal concentration of electrons 



11 T 



= 2n v , 



(A15) 



with the «,,'s denoting the subband occupations 

fh = HLk B Thih+^- B ' )/ur ]. (A16) 



When nj is fixed {i.e., the «7--constant model ), we can de- 
termine the chemical potential pi from Eq. ( |A15| ), 

7Tft 2 tlT 



m*kiiT t—k 



EM' 



+ e 



(M-S,)/k B T 



(A17) 



Appendix B: COEFFICIENTS C;'S 

Using the continuity of V g and its first derivative together 
with the (Dirichlet) boundary conditions at the end points 
Vg(-L) - V a and V g (L) - Vb, we can determine the coeffi- 
cients Cj's appearing in Eq. (46). In the regions I and V we 
find 



with 



2Ld-w V_ 

c\ wC_ - wC + , 

2L L 

C2 - ~~(2L d - w) wC- - L wC+ + V+ 

2Ld-w y_ 

Co = wC_ + wL+ , 

2L L 

cm - +^{2L d - w) wC- - LwC+ + V+ 



C ± = l -(A ± B), V ± = i(V fl ± Vt), 



(Bl) 
(B2) 
(B3) 
(B4) 

(B5) 



and A and B defined in Eq. ( A7 ). In the modulation doping 
regions 77 and IV, we have 

w 2 -2wL d + 2LL d ^ fj V. 

c-i = — C + (L d - w)C + - — 

1.. ., _ 1 



ca = +^(L d - wf C + -{L 2 d - 2wL)C + + V + , 

w 2 -2wL d + 2LL d V- 

ci = — C - (L d - w)C + - —, 

c 8 = --(Ld ~ wf C + \(L 2 d - 2wL)C + + V + . 



In the central region III, we have 
2L — 2L d + w 



cs = +- 



2L 



wC- 



c 6 = --(2L - 2L d + w) wC + + V + 



(B6) 
(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 
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